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CHAPTER  1 


INTRODUCTION 


1 . 1  Background 

Hie  addition  of  the  space  shuttle  to  the  present  booster  fleet 
will  allow  orbital  insertion  of  much  larger  structures  than  previously 
possible.  These  structures  are  generically  termed  large  space  structures 
(LSS).  The  dimensions  of  LSS  being  considered  for  future  applications 
are  on  the  order  of  20  meters  to  several  hundred  meters,  with  some 
estimates  ranging  as  high  as  several  kilometers^1 * .  Possible 
applications  range  from  telescope  systems  to  solar  power  satellites. 

LSS  are  characterized  as  being  lightweight,  flexible,  and  having 
little  natural  damping.  These  structures  tend  to  have  very  low  natural 
vibration  frequencies.  This  effect,  coupled  with  low  damping  ratios  due 
to  a  lack  of  a  viscous  medium  like  atmosphere,  causes  the  system  to  have 
long  energy  decay  times.  This  feature  makes  the  control  design  very 
difficult. 

Hie  LSS  control  problem  can  be  broken  into  three  broad 
categories;  attitude  control,  vibration  control,  and  shape  control^3!  . 

As  a  practical  matter,  simultaneous  control  of  two  or  all  three  functions 
may  be  required.  Attitude  control  for  a  satellite  or  space  vehicle  means 
the  orientation  or  pointing  of  the  rigid  body  modes,  vibration  control 
for  LSS  means  the  damping  of  vibration  modes.  Shape  control  for  LSS 
means  the  maintenance  of  a  portion  or  all  of  the  LSS  surface  contour  to  a 
prescribed  shape.  A  satellite  or  a  space  vehicle  can  be  defined  as  a  LSS 
if  the  attitude  control  system  or  other  structural  disturbances  induce 
vibrations  in  the  structure  which  are  large  enough  to  deteriorate  the 
vehicle  performance  (e.g.  line  of  sight  error)  or  in  the  more  extreme 


case  cause  shape  distortion.  Thus,  in  order  for  a  LSS  to  have  adequate 
performance  there  must  be  vibration  control. 

The  tools  of  modern  control  theory  such  as  Linear-Quadratic- 
Gaussian  (LQG)  control  theory  can  be  effectively  applied  to  the  vibration 
control  problem,  however  some  modifications  are  required^4) .  LSS  are 
essentially  distributed  parameter  systems  that  should  be  described  via 
partial  differential  equations.  In  practice,  the  dynamics  of  a  structure 
are  analyzed  using  finite  element  techniques,  this  leads  to  a  LSS  model 
which  is  composed  of  a  large  number  of  second  order  differential 
equations.  Unfortunately,  the  order  of  the  model  is  too  high  for 
practical  use  of  LQG  control  synthesis.  This  is  due  to  the  numerical 
difficulty  of  solving  the  control  equations and  the  inaccuracy  of 
the  finite  element  model  at  high  frequencies^]  .  A  common  approach  to 
overcome  this  problem  is  to  design  the  controller  using  a  reduced  order 
model,  say  by  choosing  a  set  of  modes  most  important  to  performance 
(termed  here  primary  modes) I?, 8].  Uie  controller  so  designed  is 
subsequently  applied  to  the  full  order  system.  This  leads  to  the 
undesirable  interaction  of  the  controller  with  the  ignored  modes  which  in 
turn  can  lead  to  system  instabilities^7^ .  Cie  undesirable  interaction 
of  the  controller  with  the  ignored  modes  is  often  referred  to  as 
"spillover"  effect. 

Considerable  attention  has  been  devoted  to  the  design  of  reduced 
order  controllers^9'1  .  However,  very  little  work  has  been  done  in 
the  c.rea  of  choosing  actuator/sensor  placements.  Yet  this  is  a 
fundamental  constraint  of  the  control  problem.  Essentially,  our  ability 
to  control  a  flexible  system  can  only  be  as  good  as  our  actuator/sensor 
placement  allows.  For  example,  if  actuators  and  sensors  are  located  only 
on  zero  deflection  points  of  a  mode  shape,  that  mode  cannot  be 
controlled.  This  study  undertakes  the  development  of  a  methodology  for 
actuator  and  sensor  placement  to  allow  vibrational  control  of  the  system. 


1.2 


Current  Methods  for  Actuator/Sensor  Placement 
It  is  generally  recognized  that  the  actuator  placement  problem 
and  the  sensor  placement  problem  are  dual  problems 1 J  .  This  study 
focuses  on  the  actuator  placement  issue  only.  Current  methods  for 
actuator  placement  range  from  intuitive  placement  to  quite  complex 
procedures  based  on  controllability  assessment  schemes.  Although 
intuitive  schemes  seem  to  work  quite  well  for  simple  structures  (e.g. 
beam,  plate)  they  deteriorate  as  the  complexity  of  the  structure  and  the 
order  of  the  model  increase.  For  example,  in  a  free-free  beam  modeled 
with  the  first  two  vibrational  modes  it  is  clear  that  we  want  to  locate 
actuators  near  the  peaks  of  the  mode  shapes  for  the  best  control 
authority  over  the  system.  However,  as  the  order  of  the  model  is 
increased,  it  is  unclear  which  locations  should  be  chosen  since  each  mode 
will  have  a  different  effect  on  the  system  performance  and  there  are 
multiple  peaks  in  each  mode  shape. 

There  have  been  several  attempts  to  determine  optimal  actuator 
placement  using  controllability  criteria^1  2-1 51 .  Vander  Velde  has  even 
extended  the  concept  to  include  reliability  and  redundancy  issues^15}. 

The  basic  philosophy  of  these  schemes  is  to  define  a  measure  of  the 
degree  of  controllability  and  use  it  to  evaluate  alternative  placements 
of  actuators.  The  objective  is  to  find  a  set  of  actuators  with  the 
maximum  degree  of  controllability.  Unfortunately,  these  methods  tend  to 
be  computationally  intensive  and  thus  are  not  applicable  for  high  order 
systems.  Hie  reason  for  these  computational  difficulties  is  two  fold. 
First,  the  degree  of  controllability  is  calculated  via  a  controllability 
grammian,  a  non-trivial  procedure.  Second,  as  the  number  of  possible 
actuator  locations  is  increased,  the  number  of  possible  actuator  sets 
that  must  be  evaluated  increases  dramatically.  Another  problem  with 
these  methods  is  that  they  neglect  to  take  spillover  into  account,  thus, 
a  method  for  actuator  placement  more  compatible  with  high  order  systems 
must  be  developed. 
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1 .3  Goals  of  Study 

What  we  hope  to  attain  in  the  present  effort  is  a  simplified 
method  for  actuator  placement  for  vibration  control  of  LSS.  Our  purpose 
is  to  devise  a  scheme  which  overcomes  the  computational  difficulties 
associated  with  high  order  systems.  Ideally  this  method  would  choose 
actuator  placements  such  that  they  have  good  control  authority  over  the 
primary  modes  while  talcing  spillover  considerations  into  account. 

1 .4  Approach 

It  is  possible  to  represent  both  actuators  and  linear  system 
performance  measures  as  vectors  in  modal  coordinates,  these  vectors  are 
termed  node  shapes  and  their  elements  are  the  modal  displacements 
appropriate  to  the  actuator  or  the  performance  measure.  The  actuator 
selection  problem  becomes  that  of  a  choice  of  a  set  of  actuator  node 
shapes  such  that  they  form  a  good  match  to  the  performance  node  shape. 
This  idea  can  be  stated  mathematically  as 


L,  ■  ?  •  SL 


where 

f^  is  the  performance  node  shape 
»  means  "approximately  equal" 


F  is  a  matrix  of  actuator  node  shapes 
a  is  some  arbitrary  vector  of  appropriate  dimension 


If  there  is  an  initial  set  of  N  actuator  node  shapes  (corresponding  to  N 
possible  actuator  locations)  equation  (1-1)  can  be  formulated  as  the 
following  optimization  problem 

m  m 

u  u 

min  [f  -  Y  a.f.]  Dff  -T  o.f.l  (1-2) 

-'mu  L  j-1  1  jii  3-0 J 

i.e.  find  the  set  of  mu  actuator  node  shapes  which  best  approximate 
in  the  least  squares  sense.  The  D  matrix  is  to  weight  the  modes 
according  to  their  contribution  to  performance.  To  our  knowledge  the 
above  problem  has  not  been  solved.  Thus,  two  procedures  are  suggested  to 
obtain  an  acceptable  (suboptimal)  solution  to  the  actuator  placement 
problem. 

In  the  first  method  the  initial  set  of  N  actuators  is  reduced 
using  a  QR  factorization  technique  to  solve  the  least  squares  problem. 

The  QR  factorization  as  employed  here  can  be  used  to  select  a  set  of 
independent  actuators  to  control  the  system.  In  the  second  method, 
actuators  are  selected  one  by  one,  each  actuator  complimenting  the 
already  existing  set.  Actuators  are  chosen  so  as  to  increse  the  control 
authority  of  the  actuator  set  while  minimizing  spillover. 
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1 .5  Results 

The  above  procedures  are  used  to  place  actuators  on  a  cantilever 
beam.  Both  procedures  result  in  actuator  placement  that  match  the 
intuitive  concept  of  placing  actuators  near  peaks  of  mode  shapes.  The 
effect  of  the  actuator  placement  on  the  performance  and  stability  of  an 
LQG  based  controller  is  also  demonstrated. 

1 .6  Organization  of  Thesis 

The  remainder  of  the  thesis  is  divided  into  seven  chapters. 
Chapter  2  discusses  finite  element  models  and  the  node  shape  concept. 
Chapter  3  discusses  reduced  order  controllers.  Chapter  4  presents  the 
actuator  selection  problem  and  Chapter  5  presents  two  methods  for 
actuator  selection.  Chapter  6  discusses  the  application  of  the  selection 
procedures  to  an  example  structure.  Chapter  7  contains  concluding 
remarks. 
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CHAPTER  2 


SYSTEM  REPRESENTATION 


2.1  Coordinate  Transformation 

Large  space  structures  are  lightweight,  flexible  and  lightly 
damped.  They  are  essentially  distributed  parameter  systems  that  should 
be  represented  by  partial  differential  equations.  A  convenient  model  for 
structures  is  obtained  via  an  infinite  number  of  modes  of  vibration  plus 
the  rigid  body  modes.  In  practice,  it  is  sufficient  to  analyze  these 
structures  retaining  a  finite  number  of  modes  obtained  via  the  finite 
element  method.  Computer  packages  such  as  NASTRAN  are  commonly  used  for 
this  purpose^16!. 

In  finite  element  analysis,  the  structure  is  modeled  as  a  large 
number  of  elements  connected  at  nodes1.  The  mass  of  the  structure  is 
considered  to  be  lumped  at  the  nodes  and  the  connecting  elements  are 
modeled  as  perfect  springs  and  dampers.  This  representation  is  often 
transformed  from  physical  to  modal  (or  normal)  coordinates.  This  allows 
structural  dynamics  to  be  viewed  as  a  set  of  uncoupled  modal  motions. 
Since  low  frequency  modes  will  dominate  the  motion  of  the  structure,  it 
is  reasonable  to  approximate  the  structure  by  eliminating  the  high 
frequency  modes  from  the  model. 

The  "physical"  system  equations  of  motion  obtained  via  the  finite 
element  method  with  N  nodes  is  given  by 


1 

In  finite  element  analysis  a  node  is  defined  as  a  grid  point  on  the 
structure,  not  necessarily  a  point  where  the  modal  deflection  is  zero. 
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where 


M2.  +  D3  +  *3. 


(2-1  ) 


is  a  6N  vector  of  generalized  displacements  (i.e.  each 
node  has  6  degrees  of  freedom,  3  translational  and  3 
rotational) 

is  a  vector  of  generalized  external  forces 

is  a  real  ,  symmetric,  positive  definite  mass  matrix 

is  a  real  matrix  of  inherent  structural  dancing 

is  a  real,  symmetric,  positive  semi -definite  stiffness 


The  system  equation  (2-1 )  is  conveniently  transformed  into  modal 
coordinates  via  the  transformation 


(2-2) 


where 


is  the  matrix  mode  shapes  orthonormalized  with  respect  to 
the  mass  matrix  as  in  equation  (2-4) .  (See  discussion  in 
the  following  section) . 

is  the  vector  of  modal  coordinates  (amplitudes) 


If  M  modes  are  retained  in  the  finite  element  model  then  $  is  of 
dimension  6NxM  and  is  of  length  M.  Substituting  equation  (2-2)  into 
equation  (2-1)  yields 


M$jq_  +  D«fj.  +  K$r]_  ”  h 


where  t  is  chosen  such  that 


At  this  step  D  is  defined  so  that  it  satisfies 


$t55  -  2zn 

where  Z  is  a  diagonal  matrix  of  damping  ratios  (chosen  arbitrarily) 
The  vector  of  generalized  forces  can  be  written  as 

h  =»  B  u  +  B  ,d 

—  Ur-  a— 

Thus,  the  system  equation  in  modal  coordinates  is  given  by 

— .  -2  -T-  -T- 

£  +  2Z0n  +  Jl  “  *  B^u  +  #  B^d 

where 

Z  -  DIAG(r  , ) 

i  ns 

jj  ■  OIAG(a)  ,•••,(!»  ) 

i  m 

*  *  1&,  r  •  •  *  rij 

li  is  the  mode  shape  vector  of  the  ith  mode  of  vibration 
Cl  is  the  damping  ratio  of  the  ith  mode 
u>i  is  the  natural  frequency  of  the  ith  mode 
ii  ■  (ui ,  •  •  •  the  control  inputs  from  actuators 

jj  ■  the  disturbance  inputs  from  m<j 

disturbance  sources 


r 


B  and  B  are  selection  matrices  (see  following  section)  . 


2.2  The  Modal  and  Selection  Matrices 

2.2.1  Mode  Shape 

As  stated  in  the  last  section,  the  modal  matrix  ♦  is  the  matrix 
of  mode  shapes  orthonormalized  with  respect  to  the  mass  matrix. 

Definition  1  -  'Hie  ith  mode  shape  ^  is  the  vector  of  the  displacements 
of  the  N  nodes  due  to  a  unit  amplitude  of  the  ith  mode. 

Since  there  are  six  degrees  of  freedom  for  each  node,  the  mode  shape 
vector  is  of  dimension  6N  where  the  first  N  elements  correspond  to 
displacements  in  the  first  degree  of  freedom. 

2.2.2  Selection  Matrix 

Equation  (2-2)  can  be  used  to  find  the  physical  displacements  in 
6  degrees  of  freedom  given  the  modal  coordinate  vector  r^.  Consider  the 
motion  of  a  single  node  in  one  of  its  degrees  of  freedom,  x.  Tftis  can  be 
written 

x  -  bT*Ti  (2-9) 

where  b  is  a  6N  boolean  selection  vector  with  one  non-zero  element  that 
specifies  the  node  and  degree  of  freedom  of  interest.  If  we  generalize  x 
to  a  vector  of  j  displacements,  we  can  form  the  selection  matrix  - 

B  »  [b.-] ,  •  •  • ,  b^  ]  to  obtain 


-T- 

X_  -  B  $T]_ 


(2-10) 
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2.2.3  Node  Shape 

In  analogy  to  the  concept  of  the  mode  shape,  it  is  convenient  to 

describe  the  modal  displacements  of  a  single  node  in  a  specified  degree 

-T 

of  freedom.  Thus,  define  the  product  $  5d  as  the  node  shape  vector,  f_, 
associated  with  £. 

Definition  2  -  A  node  shape  is  a  vector  of  length  M,  the  entries  of  which 
are  the  displacements  of  a  node  due  to  each  mode  excited  to  a 
unit  amplitude. 

Thus  the  nodal  displacement  in  the  specified  degree  of  freedom  can  be 
written 


x 


-T- 

The  product  B  is  defined  as  the  node  shape  matrix,  F. 

The  node  shape  is  a  versatile  concept  when  la  is  not  restricted  to 
be  a  boolean  selection  vector.  We  demonstrate  below  how  node  shapes  can 
be  used  to  define  quantities  such  as  translation  or  rotation  in  an 
arbitrary  direction,  axial  displacement  between  nodes,  and  the 
performance  of  the  system. 

2. 2. 2.1  Translation  node  shape 

Consider  the  translation  of  a  node  in  the  direction  of  a  unit 
vector  r  *  (r^,  r2»  r3, 0,0,0)  .  The  physical  translational 
displacement  vector  of  a  node  is  given  by  jc  ■  (xi ,X2»X3, 0,0,0) ,  so 
the  displacement  in  the  direction  _r  is 


x 


r 


x  •  r 


r  •  r 


x  »  r 


(2-11) 


To  form  the  dot  product  for  the  ith  node  in  modal  coordinates  we  specify 
the  selection  vector  b^  via  its  entries  b*  as  follows 


b 


k 


{r1  when  k  ■  i 

r2  k  -  N  +  i 

r3  k  -  2N  +  i 


^0  OTHERWISE 


The  resulting  node  shape  is 


It  can  be  used  to  find  the  translation  in  the  direction  r  as 


x 


r 


(2-12) 


Thus,  the  node  shape  could  be  used  to  represent  a  position  sensor  or  a 
force  actuator  oriented  in  the  direction  _r.  This  is  trivially  extended 
to  represent  a  rotation  vector  in  a  specified  direction  at  a  selected 
node. 


2. 2. 2. 2  Axial  node  shape 

The  axial  node  3hape  is  representative  f  the  axial  elongation  of 
a  member  connecting  two  nodes.  If  the  member  is  much  longer  than  the 
expected  elongation,  it  is,  approximately  correct  to  write 


x 


12 


—2 


r 


iii 


•  r 


(2-13) 


where 


12 


jc  ^and  jt 


is  the  elongation  or  contraction  of  the  element  connecting 
nodes  1  and  2 

are  the  translational  displacement  vectors  at  nodes  1  and  2 


r  is  a  unit  vector  directed  at  node  2  from  node  1 


Using  Equations  (2-11)  and  (2-12),  Equation  (2-12)  can  be  written  as 


which  shows  the  axial  node  shape,  as  given  by  the  difference 
between  the  two  translational  node  shapes  in  the  direction  of  the 
connecting  element,  i.e. 


— 1  2 


—2  ~  —1 


to  give  the  axial  elongation 


2. 2. 2. 3  Performance  node  shape 

Hie  performance  of  a  system  is  frequently  given  by  some  linear 
combination  of  the  translations  and  rotations  of  some  nodes  on  the 
structure.  A  typical  example  is  the  line  of  sight  (LOS)  error  for 
optical  systems.  Although  not  in  strict  accordance  with  definition  2,  we 
can  form  a  node  shape  that  represents  system  performance.  Considering 
the  performance  measure  y 


y  -  a,*,  +  a2x2  ♦  ...  +  Vn 


(2-14) 


where 

ai  -  an  are  constant  coefficients 

xi  -  xn  are  translational  or  rotational  deflection  of  nodes 
Equation  (2-14)  can  be  rewritten  as 
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CHAPTER  3 


LARGE  SPACE  STRUCTURE  CONTROLLERS 


3.1  Introduction 

As  noted  in  Chapter  2,  finite  element  methods  can  be  used  to 
approximate  a  system  of  infinite  order  using  only  M  modes.  Unfortunately 
the  dimension  of  the  system  is  still  too  large  to  allow  construction  of  a 
practical  full  state  feedback  controller.  This  is  due  to  onboard 
computer  limitations  and  the  difficulty  in  computing  optimal  control 
gains  for  high  order  systems^],  a  common  approach  to  the  control 
design  is  to  identify  a  set  of  primary  modes,  usually  based  on  system 
performance,  and  design  a  controller  based  on  those  modes  only^7'8!. 

The  goal  is  the  control  of  the  full  order  system  with  this  reduced  order 
controller. 

4 

3.2  State  Equation  Development 

The  system  equation  is  given  in  Equation  (2-8)  as 

..  — •  -2  -T-  -T- 

r±  +  2ZJ2_q^  +  (5  *»  4  B^u  +■  4  B^^d 

In  addition,  the  measurement  equations  can  be  written  as 

**  -  6x*a  +  ^  t3-1) 


21 


[i 


*V  "  V*  +  ^ 


(3-2) 


where 


y*  is  a  vector  of  position  sensor  output 
is  a  vector  of  velocity  sensor  output 
and  are  the  appropriate  selection  matrices 

Vjj  and  Vy  represent  sensor  noise 


The  corresponding  state  equation  is 


Ax  +  Bu  +  Td 


2_  m  Csc  +  v 


(3-4) 


For  the  representation  in  (2-8),  (3-1)  and  (3-2),  the  various 
entries  of  (3-3)  and  (3-4)  are  defined  as  follows: 


-2  *  — 
-U  I-2ZJ 


v 

—x 

V 

—V 

.  4 


-T 

B  $  1  0 

_*..l _ 

I  -x 

0  I  B  « 

■  tr  1 


By  selecting  L  primary  modes,  equations  (2-8),  (3-1),  and  (3-2) 
can  be  represented  via  primary  and  residual  components. 

‘ip  +  2Vp%  +  ~  *p5up^  +  *p5dpd 

Hr  +  ^rVr  +  4lR  “  ^  +  iR5d*£ 

^X  *  VpD,  +  ®xR*R^R  +  ~X 
■^v  “  \pVp  +  ®vR®R^R  +  ~v 


where  subscript  P  is  associated  with  primary  modes  and  subscript  R  with 
residual  modes. 

Equations  (3-6)  are  written  in  state  equation  form  as 


(3-7) 


3.3  Controller  Design 

The  reduced  order  controller  is  designed  based  on  the  primary 
modes  only.  Restricting  attention  to  optimal  (modern)  control 
philosophy,  the  controller  is  comprised  of  two  parts: 

1)  A  state  estimator  which  accepts  the  sensor  measurements  £ 

A 

and  produces  an  estiamte  x  of  the  primary  states  x 


2)  A  linear  state  feedback  control  law  which  multiplies  the 

A 

state  estimate  by  a  constant  gain  matrix  to  produce  the 
actuator  commands  u. 
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The  state  estimator  equation  is 


x„ 
— P 


A  >c 

P-P 


+  BpU  +  G  -  0^) 


(3-8) 


where  G  is  the  constant  estimator  gain  matrix.  The  control  law  is  given 
by 


u  =* 


(3-9) 


where  K  is  the  constant  feedback  gain  matrix.  Using  equations  (3-7)  and 
(3-9)  we  can  rewrite  (3-8)  as 


^P 


(Sp 


V 


-  5cn]x 

pJ—p 


+  GC  x 
P—P 


+  50  * 
R-R 


+  Gv 


(3-10) 


Defining  the  estimate  error,  e_  as 


e  =  x„  -  x„ 
—  — P  — P 


The  estimate  error  equation  can  be  written  as 


(Ap  -  GCph 


GC  x 
R— R 


So  the  total  system  state  equations  can  be  written  as 
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Block  diagrams  of  the  controller  and  the  controlled  system  are  given  in 
Figures  (3-1)  and  (3-2). 


Figure  3-1  Feedback  Controller 


The  selection  of  the  gain  matrices  K  and  G  is  a  standard 
procedure  given  the  correct  plant  model^1^.  However,  the  task  is 
complicated  due  to  the  reduced  order  model  used  in  the  controller 
design,  which  in  turn  leads  to  spillover  effect  (discussed  in  next 
section) .  The  proper  selection  of  gain  matrices  to  accommodate  spillover 
is  beyond  the  scope  of  this  study  and  will  not  be  discussed  (see  [9-10] 
for  some  possible  solutions) . 
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3.4  Spillover 

One  of  the  major  problems  to  be  considered  in  reduced  order 
controller  design  is  spillover.  Spillover  is  the  adverse  interaction  of 
the  controller  with  the  residual  modes.  Considering  Figure  3-2,  the 
residual  modes  are  excited  by 


B  u  =  B  Kx. 
ft—  R  — P 


FORCE  ACTUATORS 


FLEXIBLE  SYSTEM 


SENSORS 


MEASUREMENT 

NOISE 


Figure  3-2  System  with  Reduced  Order  Controller 

This  is  called  control  spillover.  Furthermore,  the  measurements,  y  are 
contaminated  by  the  term  C _x_.  'Biis  in  turn  lead  to  errors  in  the 
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m 


estimate  The  contamination  of  sensor  readings  by  residual  modes  is 

termed  observation  spillover. 

Control  spillover  alone  is  probably  not  a  major  problem.  If  there 

is  no  observation  spillover  (i.e.  C  -0),  then  from  equations  (3-11)  the 

R 

eigenvalues  of  the  controlled  system  can  be  obtained  from  (Ap+BpK), 

(Ap-GCp),  and  Ap.  Assuming  the  residual  modes  are  open  loop  stable 
(even  if  only  marginally  stable)  and  that  the  regulator  and  observer  are 
stable,  the  only  effect  of  control  spillover  is  to  excite  the  residual 
modes.  This  will  tend  to  deteriorate  performance,  but  will  not  lead  to 
instability. 

However,  observation  spillover  might  lead  to  instabiltyf7 J . 

Considering  Equation  (3-11),  if  C  *  0  then  the  system  eigenvalues  will 

R 

be  dependent  on  C  .  Thus  the  closed  loop  poles  might  migrate  to  the 
right  half  s-plane,  especially  considering  the  low  damping  of  the 
structure. 
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CHAPTER  4 


ACTUATOR  PLACEMENT  PROBLEM 


4.1  Introduction 

Actuator  locations  are  a  fundamental  constraint  of  the  LSS 
controller  design  problem.  It  is  desirable  to  place  actuators  such  that: 

1 )  control  authority  over  the  L  primary  modes  is  maximized 

2)  spillover  to  residual  modes  is  minimized 

3)  a  small  number  of  actuators  is  used 

4)  the  actuator  set  is  balanced,  namely  that  actuators  are 
operating  at  comparable/acceptable  power  levels 

The  actuator  placement  problem  is  usually  subject  to  the  following 
constraints 


1 )  actuator  type 

2)  physical  limitations  on  actuator  location 

In  this  chapter  the  constraints  and  placement  criteria  are  discussed  in 
detail,  then  a  mathematical  formulation  of  the  actuator  placement  problem 
is  presented. 


4.2  Constraints 


4.2.1  Actuator  Type 

In  most  design  applications,  the  choice  of  actuator  type  is 
constrained.  This  might  be  due  to  problems  with  contamination  by 
propellants,  power  supply  limitations,  cost,  etc.  An  important 
consideration  in  selecting  actuators  for  LSS  applications  is  the 
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availability  of  broad  band  (especially  low  frequency)  actuators.  Another 
desirable  actuator  characteristic  is  its  operation  with  no  interaction 
with  the  rigid  body  modes.  Thus  propellant  devices  and  momentum  wheels 
are  of  limited  use  in  vibration  control.  Member  actuators  are  desirable 
due  to  their  lack  of  effect  on  the  rigid  body  motion  and  their  spectral 
characteristics. 

4.2.2  Physical  Limitations 

There  are  locations  on  a  structure  where  the  placement  of 
actuators  is  forbidden.  These  forbidden  regions  are  usually  due  to 
1)  the  sensitivity  of  nearby  equipment  or  2)  portions  of  the  structure 
being  too  weak  to  withstand  actuator  inputs  of  the  magnitude  required. 

In  the  first  case,  effects  like  the  outgassing  of  propellant  actuators 
must  be  considered.  In  the  second  case,  locations  such  as  mirror  mounts 
are  usually  forbidden. 

4.3  Criteria  for  Placement 

4.3.1  Control  Authority 

Actuator  placement  to  maximize  the  control  authority  over  the 
primary  modes  means  selecting  actuator  locations  such  that  the  actuators 
are  operating  at  the  minimum  possible  power  levels  to  meet  performance 
specifications.  Por  example,  consider  a  single  mode  shape.  It  is  clear 
that  locating  an  actuator  at  the  peak  of  the  mode  shape  will  give  better 
control  authority  over  that  mode  than  an  actuator  located  off  the  peak. 
Thus,  actuator  node  shapes  that  provide  good  control  authority  should 
have  large  elements  corresponding  to  the  primary  modes. 

4.3.2  Spillover 

It  is  desirable  to  place  actuators  such  that  excitation  of  the 
residual  modes  is  minimal.  Wiis  in  turn  would  result  in  better  control 
when  the  reduced  order  controller  is  applied  to  the  structure.  There  are 


tvw  considerations  when  choosing  actuator  locations  to  minimize 
spillover.  First,  actuators  should  be  located  as  close  to  the  zero 
deflection  points  of  the  residual  mode  shapes  as  possible.  In  general, 
this  is  not  possible  for  all  of  the  residual  modes.  A  possible  solution 
is  to  choose  a  subset  of  the  most  critical  residual  modes,  termed 
secondary  modes,  and  choose  actuator  locations  near  their  zero  deflection 
points  (i.e.  small  entries  in  the  actuator  node  shapes  corresponding  to 
residual  modes).  The  second  consideration  is  the  applied  actuation 
force.  If  an  actuator  is  moved  away  from  a  point  where  it  has  good 
control  authority,  it  will  require  larger  actuation  forces  to  control  the 
system.  In  general,  larger  actuation  forces  imply  more  spillover  to  the 
residual  modes.  Thus,  there  is  a  tradeoff  to  be  considered  in  choosing 
actuator  locations  for  spillover  prevention. 

4.3.3  Number  of  Actuators 

Choosing  the  optimum  number  of  actuators  to  control  the  system  is 
not  a  trivial  problem.  For  economic  considerations,  the  number  of 
actuators  should  be  kept  small.  However,  if  redundancy  is  a  major 
consideration,  the  actuator  set  should  be  larger  than  the  minimum 
required  to  maintain  control.  For  the  purposes  of  this  study,  redundancy 
issues  are  ignored  and  the  smallest  set  that  retains  good  control 
authority  is  sought. 

A  system  is  said  to  be  controllable  if  it  can  be  transferred  from 
any  initial  state  x.(0)  to  a  final  state  x,(tf)  in  finite  time  by  some 

control  n.  Thus,  controllability  is  the  ability  of  the  control  input  to 

affect  each  state  variable.  Now,  consider  a  LSS  with  L  primary  modes. 

The  above  implies  that  a  single  actuator  can  control  the  primary  modes  if 
all  the  corresponding  elements  in  the  actuator  node  shape  are  non-zero 
(note,  this  is  true  as  long  as  each  mode  has  different  frequencies) .  On 
the  other  hand,  L  actuators  can  be  found  such  that  practically  each  mode 
is  independently  controlled  by  a  single  actuator. 
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Unfortunately  controllability  is  simply  a  binary  concept  -  either 
a  system  is  controlled  or  it  is  not.  It  is  intuitively  obvious  that 
increasing  the  number  of  actuators  will  increase  the  control  authority  of 
the  actuator  set,  just  as  increasing  the  number  of  measurements  will 
increase  the  accuracy  of  an  estimte  produced  by  an  estimator.  However, 
if  initial  actuator  locations  are  well  selected,  diminishing  returns  are 
expected  as  the  number  of  actuators  approaches  L^l  .  if  actuators  are 
placed  such  that  individual  actuators  have  good  control  authority  over 
several  modes,  then  it  is  possible  that  less  than  L  actuators  will 
provide  approximately  the  same  control  authority  over  the  primary  modes 
as  L  actuators. 

4.3.4  Balanced  Set 

It  is  desirable  to  have  the  actuators  working  at  comparable  power 
levels.  This  is  convenient  since  identical  actuators  can  be  used 
throughout  the  structure  which  simplifies  the  design  procedure  and 
reduces  cost. 

4.4  Mathematical  Formulation 

Expressing  the  above  objectives  in  a  mathematical  expression 
which  is  ameanable  to  actuator  placement  procedures  is  difficult  if  not 
impossible.  In  the  following,  some  possible  engineering  solutions  which 
approach  the  goals  stated  above  are  presented. 

In  Chapter  2,  the  concept  of  the  node  shape  is  developed.  It  is 
demonstrated  that  node  shapes  can  be  created  to  represent  actuator 
influence  and  system  performance  measures.  Thus,  the  problem  of  actuator 
placement  can  be  mathematically  stated  as  that  of  choosing  a  set  of 
actuator  node  shapes  such  that  a  linear  combination  of  them  is 
approximately  equal  to  the  performance  node  shape  i.e. 

-  F  •  a  (4-1) 
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where 

_f*  is  the  performance  node  shape 

-  means  "approximately  equal" 

F  is  a  matrix  of  actuator  node  shapes 

a  is  some  arbitrary  vector  of  approximate  dimension 

The  above  formulation  can  be  somewhat  simplified  if  an  initial 
set  of  n  possible  actuator  locations  has  already  been  selected.  This 
initial  set  is  selected  subject  to  constraints  on  actuator  type  and 
location.  The  problem  then  reduces  to  a  selection  of  my  node  shapes 
from  the  initial  set  which  provide  a  good  approximation  of  f_* .  This 
selection  is  subject  to  the  criteria: 

1 )  Maximize  control  authority:  The  match  to  the  performance 
node  shape  is  good,  entries  in  F  corresponding  to  primary 
modes  are  large,  <x  components  are  small 

2)  Minimize  spillover:  The  entries  in  F  corresponding  to 
secondary  modes  are  small. 

3)  Small  actuator  set:  my  is  small 

4)  Balance:  All  the  a:  components  are  of  comparable  magnitude 

In  the  following  section,  the  quadratic  cost  function  is  proposed 
as  the  mathematical  criterion  for  "approximately  equal". 

4.4.1  Lea3t  Squares  Formulation 

The  selection  criteria  can  be  formulated  as 

m 

|i|2  •  Ll  -  l  u  a/  Hi. 

11  i-1 


l  £&) 

i-1 


(4-2) 
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where  mu  and  a  are  chosen  to  minimize  the  squared  magnitude  of  the 
residual  vector  r^.  Hie  purpose  of  the  diagonal  matrix  5  is  to  weight  the 
different  modes  according  to  their  importance  to  the  control  oblem.  In 
effect,  the  5  matrix  reduces  the  order  of  the  problem  to  conL^der  only 
the  primary  modes.  Thus,  L  non-zero  values  and  m-L  zero  values  are 
expected  among  the  elements  of  D.  A  possible  selection  for  D  is  to  use 
ones  for  the  elements  corresponding  to  residual  modes.  However,  this 
does  not  reflect  relative  importance  among  the  primary  modes.  Thus,  the 
elements  of  5  corresponding  to  primary  modes  should  reflect  the 
contribution  of  each  mode  to  overall  performance. 

For  example,  for  a  given  disturbance  acting  on  the  system,  a  set 

of  mean  square  modal  amplitudes,  rima >  can  be  obtained.  If  the  intra 

mode  interaction  is  small,  it  is  reasonable  to  approximate  the  contri- 

2 

bution  of  the  ith  mode  to  the  mean  square  performance  as  (f#).(n  ). • 


5  »  diag( (f#)?)  diag  (<hms)i)  diag  ( w i> 


where 

i  ranges  from  1  to  M 
w-jL  *  1  for  primary  modes 
^  *  0  for  residual  modes 


CHAPTER  5 


ACTUATOR  SELECTION  PROCEDURES 


5.1  General  Method 

In  Chapter  4  a  set  of  constraints  and  criteria  were  developed  for 
the  actuator  placement  problem.  The  placement  problem  was  then 
mathematically  formulated  as  a  least  squares  problem.  To  our  knowledge, 
the  least  squares  problem  as  presented  has  not  been  solved,  thus  a 
suboptimal  approach  is  taken  to  solve  the  least  squares  problem.  The 
suggested  general  procedure  consists  of  the  following  major  steps: 

1 )  Define  a  large  initial  set 

2)  Choose  L  actuators  to  solve  the  least  squares  problem 

3)  Reduce  the  set  of  L  actuators  if  possible 

4)  Augment  the  initial  set  if  necessary 

rf  flow  chart  of  the  general  procedure  is  given  in  Figure  5-1 .  Each  of 
the  above  steps  is  discussed  in  detail  in  the  following  sections. 

5.2  Initial  Set 

The  initial  set  of  actuator  locations  is  ideally  an  infinite  set 
containing  all  permissable  actuator  locations.  This  obviously  being 
impossible,  a  large  set  that  is  representative  of  the  infinite  set  is 
selected.  The  first  constraint  that  must  be  considered  is  what  type  (or 
types)  of  actuators  can  be  used  to  control  the  structure.  As  noted  in 
Chapter  4,  it  is  desirable  to  use  member  actuators  since  they  do  not 
affect  rigid  body  modes.  Thus,  whenever  possible,  member  actuators 
should  be  used  in  the  formulation  of  the  initial  set.  Actuators  of  the 
type(s)  chosen  are  placed  at  a  large  number  of  locations  on  the 
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ACTUATORS 


structure,  avoiding  those  areas  where  actuator  placement  is  physically 
prohibi ted. 

For  example,  if  only  axial  actuator  are  to  be  used  to  control  a 
truss  structure,  a  simple  initial  set  would  be  the  placement  of  an 
actuator  on  each  member. 

5.3  Least  Squares  Solution  Based  on  L  Actuators 

Hie  solution  of  the  least  squares  problem  presented  in  Chapter  4 
involves  the  selection  of  my  actuators  and  the  components  of  £  to 
minimize  the  norm  of  the  least  squares  residual  vector.  As  a  first  step 
in  the  suboptimal  solution,  L  actuators  are  chosen  from  the  initial  set 
and  the  components  of  calculated.  In  this  section  a  numerical 
technique  to  solve  for  the  _a  components  and  two  methods  for  selecting  L 
actuators  from  the  initial  set  are  presented. 

Before  continuing  with  this  section,  some  new  notation  is 
required  that  willbe  used  throughout  the  chapter.  The  least  squares 
formulation  presented  in  Chapter  4  was  reduced  to  consideration  of 
weighted  primary  modes.  Thus,  the  quantities  of  interest  can  be 
represented  as 

-Dy  f 

(5)1  f 

(5) 2  i. 
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where  it  is  assumed  that  the  primary  modes  have  been  pivoted  to  the  top 
of  the  appropriate  vectors.  Thus  Equation  (4-2)  is  rewritten 


[f.  -  y  f  . a. ]T  [f, 

L— *W  .  L  — Wl  L— *w 

1-1 


5.3.1  £R  Factorization 

The  QR  factorization  is  a  technique  commonly  used  to  obtain  the  a_ 

components  in  a  least  squares  problem.  Ihe  concept  is  presented  in  the 

context  of  node  shapes,  but  it  can  be  generalized  to  any  set  of  vectors. 

Let  F„be  an  L  X  m  node  shape  matrix.  Assume  that  L  >  m  and  rank  (F  )  = 
w  u  —  u  w 

-  -T- 

m  .  Wiere  exists  an  orthogonal  matrix  Q  such  that  Q  F  is  zero  below  its 

U 

-T- 

diagonal.  Thus,  Q  pwcan  be  written  in  the  form 

(5-1 ) 


where  R  is  a  m  xm  upper  triangular  matrix.  In  the  case  when  m 
u  u  u 

Equation  (5-1 )  is  written 


L, 


-T- 

Q  F  -  R 

*  w 

which  implies 

Fw  -  QR  (5-2) 
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However,  in  the  case  where  m^  <  L,  Q  can  be  partitioned  as 


Q  =  [Qr  Q2] 


where  Q.  is  L  x  m  .  Thus,  F  can  be  written  as 
1  u 


F  =  Q  R 
w  1 


(5-3) 


Now,  the  least  squares  problem  is 


F  a  =  f* 
w —  — *w 


(5-4) 


which  can  be  rewritten  as 


V*  *  -*w 


(5-5) 


If  Fwhas  linearly  independent  columns,  as  initially  assumed,  then  R  is 


nonsingular  and 


— 1  -T  , 

2  m  R  Q1  i*w 


(5-6) 


It  can  also  be  shown  that  the  least  squares  residual  vector,  r _w,  is 
given  by 


— w  “  ®2®2^-*w 


*TI'~**-  r~ 


Thus,  in  the 


(i.e.  the  projection  of  _f ^  orthogonal  to  the  span  of  F  ) . 
case  mu=L,  j rw |  *  0. 

5.3.2  Maximum  Norm  Method 

'flie  QR  factorization  and  least  squares  routines  used  in  this 

study  are  standard  LINPACK  routines^!  .  a  pivoting  option  is  included 

in  the  QR  factorization  routine  that  appears  useful  as  a  method  to  reduce 

the  initial  set  to  L  actuators.  Given  an  actuator  node  shape  matrix  for 

the  initial  set,  F  _  of  dimension  Lx n  where  n  is  the  number  of  actuators 
w  X 

in  the  initial  set,  the  pivoting  routine  selects  a  set  of  independent 
actuators  from  the  initial  set  and  pivots  them  to  the  front  of  F 

wl 

Conceptually,  the  pivoting  procedure  is  as  follows: 

(1)  Given  FwI  =  [  f^  ,  •••,  fj,  compute  |  (fj  | 
for  all  n  actuators  (node  shapes) . 


(2) 


Pivot  node  shape  with  maximum  f  | 


to  first  column  F  „ 
wl 


(3)  Let  i  represent  number  of  actuators  already  selected.  For 
the  n-i  remaining  node  shapes,  calculate  the  component 
orthogonal  to  the  span  of  the  i  selected  node  shapes,  f-L 


(4) 

(5) 


for  the  n-i  node  shapes 


Pivot  actuator  with  maximum 

column  of  F  , 
wl 


to  the  (i+1 ) th 


(6)  If  L  node  shapes  have  been  selected,  stop,  otherwise 

return  to  step  3  and  begin  process  to  select  another  node 
shape. 
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The  QR  factorization  procedure  can  now  be  performed  on  the  first  L 
columns  of  F  .  In  practice,  the  pivoting  and  computation  of  the 

[21] 

elements  of  Q  and  R  takes  place  simultaneously  .  They  are  separated 
here  to  make  the  presentation  of  the  concepts  clearer. 

As  a  result  of  the  pivoting  option,  the  R  matrix  has  diagonal 
elements,  r^,  which  obey 


11  - 


r22|  ^ 


-  Tll 


where  |rii|  is  t^le  norm  of  the  component  of  orthogonal  to  the 
span  of  [jWl ,  fW2,  •«*,  fWi_i ] .  If  at  some  point  | r j  =0 , 
there  is  not  a  set  of  L  independent  actuators  in  the  initial  set. 

This  pivoting  routine  appears  useful  for  actuator  selection  for 
several  reasons.  First,  maximum  norm  values  are  selected,  which 
corresponds  to  the  concept  of  attaining  good  control  authority.  Second, 
on  independent  set  is  selected,  which  corresponds  to  the  concept  of  a 
small  actuator  set.  Spillover  and  balance  criteria  have  not  been 
considered  with  this  procedure. 

5.3.3  Residual  Matching  Method 

The  second  method  for  selecting  L  actuators  from  the  initial  set 
is  conceptually  similar  to  the  maximum  norm  method.  It  is  formulated  to 
select  an  independent  set  of  actuators  that  have  good  control  authority 
over  the  primary  modes  and  avoid  spillover  into  selected  secondary 
modes.  This  method  is  presented  in  the  following  two  sections:  in  the 
first  section  control  authority  is  considered  and  in  the  second  section 
spillover  is  added  to  the  formulation. 


5. 3. 3.1  Control  authority  considerations 

A  diagram  of  the  two  vector  least  squares  problem  is  given  in 
Figure  5-2. 
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To  obtain  good  control  over  the  system  with  a  single  actuator,  the  angle 
0  between  the  performance  node  shape,  and  the  actuator  node  shape, 

f_w,  should  be  small.  In  addition,  to  control  with  small  actuation 
forces,  the  length  of  £w  should  be  large.  Thus,  a  good  measure  of  the 
control  authority  of  a  single  actuator  would  be  the  length  of  the 
projection  cf  the  actuator  node  shape  onto  the  performance  node  shape. 
This  projection  length  can  be  written 


— *w  *  — *w 


Define  a  control  authority  ranking  measure,  c^,  for  the  ith  actuator  as 


the  above  definition  is  adequate  for  non-vibrating  systems,  but 
must  be  modified  for  vibrational  systems.  The  performance  of  vibratory 
systems  is  often  measured  in  terms  of  mean  square,  or  RMS,  quantities. 
In  node  shape  notation,  the  system  performance  is  measured  by 


y 


The  appropriate  mean  square  measure  is  obtained  from 

B(y2]  -  B(Uln)2]  -  E[(V,  *  —  ♦  V„)2] 


If  the  intra-mode  interaction  is  assumed  small,  it  is  approximately 
correct  to  write 


+  f. 


EW 


The  purpose  of  the  above  is  to  demonstrate  that  the  mean  square 
performance  is  not  significantly  affected  by  the  signs  of  the  elements  of 
f*  (i.e.  they  are  always  squared  quantities).  However,  the  control 
authority  ranking  measure  developed  above  is  highly  dependent  on  the 
signs  of  the  terms  in  f_*.  For  example,  if  f_ and  f_w  are  given  by 


The  ranking  scalar,  c^,  is  equal  to  zero.  This  implies  that  £w  has 
no  control  authority  over  _f*v.  This  is  only  true  in  the  case  where  the 
two  modes  have  the  same  frequency  and  are  in  phase.  In  general,  this 
will  not  be  true.  Thus,  the  control  authority  ranking  scalar  must  be 


k2 


redefined  so  as  to  ignore  the  signs  of  the  components  of  f^.  The 
ranking  scalar  is  redefined  such  that  the  absolute  value  of  each  term  in 
the  dot  product  is  used.  Therefore 


c. 

1 


L 


I 

j-1 


f  f*. 
w ,  *w . 


With  the  above  formulation,  c^  should  be  a  reasonable 
indication  of  the  control  authority  of  the  ith  actuator.  Selection  of 
the  actuator  with  the  largest  c^  value  would  be  a  reasonable  choice  as 
the  best  actuator  to  control  the  system. 

After  selecting  the  first  actuator  to  control  the  system,  it  is 
clear  that  if  j_r_w |  *  0  the  match  to  the  performance  node  shape  can  be 
improved  by  adding  an  actuator  node  shape  in  the  direction  r^,.  If  the 
remaining  actuators  in  the  initial  set  are  ranked  by 


the  best  actuator  in  the  direction  _rw  can  be  selected.  Each  time  an 
actuator  is  added  to  the  set  a  new  residual  vector  is  calculated.  Hie 
selection  process  continues  until  L  actuators  are  selected.  A  flow  chart 
of  the  procedure  is  given  in  Figure  5-3. 

A  property  of  the  residual  vector  is  that  it  is  orthogonal  to  the 
space  spanned  by  the  actuator  node  shapes.  Thus,  an  independent  set  of 
actuators  is  expected  from  this  selection  process.  Thus  the  selected  set 
should  be  an  independent  set  with  good  control  authority  characteristics. 


5. 3. 3. 2  Spillover  considerations 

The  residual  mode  equation  present  in  Chapter  3  is 
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Figure  5-3  Residual  Matching  Procedure  (Control  Authority  Only) 


^R  +  2ZRnRnR 


-T-  -T- 

<t_B  „u  +  i  B,  d 
R  uR—  R  dR— 


The  best  spillover  characteristics  would  be  obtained  if 


0 


In  general,  the  global  prevention  of  spillover  to  residual  modes  via 

actuator  placement  is  not  possible.  Thus,  a  reduced  set  of  ng  modes  is 

defined  and  actuator  locations  are  selected  to  reduce  spillover  to  those 

modes.  The  modes  in  this  set  are  termed  secondary  modes.  If  B  is 

us 

defined  as  the  rows  of  corresponding  to  secondary  modes,  then  the 
actuator  placement  problem  is  defined  in  terms  of  minimizing  B  g . 

Each  column  of  5  is  the  secondary  mode  elements,  s.,  of  an 
us  1  y 

actuator  node  shape.  Define  a  spillover  ranking  scalar,  d^,  for  the 
ith  actuator  in  the  initial  set  as 


Thus,  actuators  with  the  smallest  secondary  mode  elements  can  be  chosen 
from  the  initial  set.  This  spillover  criteria  can  be  combined  with  the 
control  authority  criteria  developed  above. 

The  suggested  procedure  is  to  first  sort  the  initial  set  by 
control  authority  characteristics  and  then  select  a  small  set  of  the 
actuators  with  the  best  control  authority.  This  set  is  then  sorted  by 
the  spillover  criteria  and  a  single  actuator  chosen.  In  this  manner  a 
set  of  independent  actuators  with  good  control  authority  can  be  obtained 
while  taking  spillover  into  consideration. 
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The  selection  of  the  small  set  from  the  initial  control  authority 
ranking  is  as  follows.  Define  c^x  as  the  maximum  control  authority 
scalar  obtained  in  the  ranking  procedure.  If  1  is  defined  as  a  limiting 
factor,  then  those  actuators  with  c^  values  less  than  lc^x  can  be 
eliminated  from  the  selection  process.  This  reduced  set  can  now  be 
sorted  by  spillover  characteristics.  In  working  with  several  example 
structures  (see  Chapter  6),  a  value  of  1=0.9  appears  to  be  a  goodchoice. 

A  flow  chart  of  the  modified  selection  procedure  is  given  in  Figure  5-4. 


5.4  Reducing  the  Number  of  Actuators 

A  desirable  set  of  actuators  achieves  good  control  authority  and 
minimum  spillover  with  a  small  set  of  actuators.  In  the  above,  a  set  of 
L  actuator  node  shapes  is  chosen  which  results  in  a  perfect  match  to  the 
performance  node  shape  (i.e.  | r_w j  =  0).  In  this  section  a  method  is 
presented  to  identify  which  of  the  L  actuators  have  the  smallest 
contribution  to  the  least  squares  fit.  It  is  possible  that  these 
actuators  can  be  removed  from  the  set  without  causing  a  significant 
increase  in  J r^w |  (i.e.  |r.w|  ■  0). 

Equation  5-5  can  be  rewritten  as 


a  =  Ra  =  oj f,w 

Hie  columns  of  =  [c^ , •  •  •  'finy  ]  form  an  orthonormal  basis  for  the 

apace  spanned  by  F^.  Note  that  g^  ■  f q^f ^  j  is  the  length  of  the 

projection  of  f _*w  on  .  Thus,  the  magnitude  of  gj.  is  an 
indication  of  the  importance  of  the  ith  actuator  to  the  least  squares 
approximation.  Hie  actuator  elimination  procedure  is  as  follows: 

1)  rank  actuators  by  g^  values 

2)  eliminate  actuator  with  lowest  g^  value 

3)  calculate  the  least  squares  residual  vector,  jrw,  for 
the  remaining  set 
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Figure  5-4  Residual  Matching  Procedure  With  Spillover  Incorporated 


4)  if  I rw I  is  less  than  some  e  then  eliminate  another 
actua tor 

5)  if  |rw|  is  greater  than  some  e  then  add  last  actuator 
eliminated  back  to  set  and  stop. 

e  is  the  criteria  for  |_r_w J  ■  0.  In  this  study  e  is  chosen  as 
1.0x10“®,  this  being  the  approximate  accuracy  of  single  precision 
computations. 

5.5  Augmenting  the  Initial  Set 

In  the  actuator  reduction  procedure  above,  it  is  possible  to 
identify  modes  which  are  poorly  controlled.  Inspection  of  the  terms  in 
the  residual  vector,  ry,  will  indicate  poorly  controlled  modes.  If 
some  of  the  terms  in  rw  are  greater  than  e,  then  no  reduction  of  the 
set  will  be  possible.  Thus,  it  is  desirable  to  find  actuator  locations 
that  will  provide  adequate  control  of  these  modes.  These  locations  are 
added  to  the  initial  set  and  the  selection  procedure  reaccomplished. 

The  selection  of  the  locations  for  such  actuators  can  be  aided  by 
a  visual  inspection  of  the  vibrating  structure  simulation  (if  one  is 
available)  and/or  the  analysis  of  the  mode  shape(s).  For  example, 
locatins  at  the  peak  of  a  mode  shape  should  provide  good  control 
authority  over  that  mode. 

5.6  Multiple  Performance  Criteria 

In  many  cases  two  or  more  linear  performance  criteria  must  be 
simultaneously  met.  For  example,  in  an  optical  structure,  the 
performance  specification  may  be  to  minimize  line  of  sight  (LOS)  error  in 
two  directions,  say  x  and  y.  This  situation  necessitates  some  minor 
modifications  to  the  previous  development. 

First,  the  primary  modes  must  be  defined  in  terms  of  all  the 
performance  criteria.  This  leads  to  a  reevaluation  of  the  weighting 
matrix  5.  Recall  that  the  elements  of  the  diagonal  D  matrix  were  the 
approximate  contribution  of  each  mode  to  the  overall  performance.  Thus, 
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an  overall  performance  criteria  must  be  defined.  In  the  LOS  error 
example  above,  overall  LOS  error  can  be  written 


LOS 


(L0Sy)2 


If  the  intra-mode  interaction  is  small,  it  is  reasonable  to  approximate 
the  contribution  of  the  ith  mode  to  the  mean  square  performance,  1^  as 


1. 

1 


U*x)J  kJi  +  u.,)?  (Oi 


where  JE*X  and  f^y  are  the  performance  (LOS)  node  shapes  for  the  x  and 
y  directions  and  nmg  is  a  vector  of  mean  square  modal  amplitudes  for 
the  uncontrolled  system.  Thus  the  elements  of  the  5  matrix,  d^,  are 
given  by 


ii 


1 .  w. 
l  i 


where  w^*1  for  primary  modes  and  0  for  residual  modes. 

No  modifications  are  required  to  use  the  maximum  norm  selection 
procedure.  However,  several  modifications  are  required  to  use  the 
residual  matching  procedure.  The  residual  matching  procedure  is  designed 
to  select  L  actuator  node  shapes  to  match  each  performance  node  shape. 
Thus,  as  many  as  2L  actuators  could  be  selected.  Some  duplication  in  the 
actuators  selected  to  ma^tch  each  performance  measure  is  expected  and  in 
some  cases  the  selections  would  be  identical.  In  those  cases  where  more 
than  L  actuators  are  selected,  dependent  actuators  must  be  eliminated 
from  the  set.  The  suggested  procedure  is  to  reduce  this  set  via  the 
maximum  norm  technique. 


Hie  final  modification  is  in  the  actuator  set  reduction 
technique.  Here,  the  elements  of  g  must  be  evaluated  for  each 
performance  criteria  before  an  actuator  can  be  chosen  for  elimination. 


CHAPTER  6 


100  METER  SCE  BEAM  EXAMPLE 


6.1  Physical  Description 

The  first  actuator  placement  study  is  devoted  to  a  100  meter 
space  construction  experiment  (SCE)  beanJ1^  shown  in  Figure  6-1.  Hie 
SCE  beam  is  designed  to  be  deployed  from  the  space  shuttle  cargo  bay  and 
will  be  a  test  bed  for  various  LSS  experiments.  The  SCE  beam  is  of 
diamond  truss  construction,  is  rigidly  attached  to  the  space  shuttle,  and 
has  a  250  kg  tip  mass.  Additional  properties  are  given  in  Appendix  A. 

6.2  Finite  Element  Analysis 

Hie  beam  is  modeled  as  having  61  elements  connected  at  54  node 
points.  Node  locations  and  element  connections  data  are  tabulated  in 
Appendix  A.  Hie  finite  element  model  is  based  on  10  vibrational  modes 
and  6  rigid  body  modes.  Open  loop  vibration  frequencies  are  listed  in 
Table  6-1  and  plots  of  the  mode  shapes  are  contained  in  Appendix  A. 

6.3  Disturbance  Description 

For  the  purposes  of  this  study,  two  white  noise  disturbances  are 
applied  to  node  39,  one  is  oriented  in  the  x  direction  and  the  other  in 
the  y  direction.  Each  disturbance  is  zero  mean  and  has  a  variance  of 
800  nt2. 


6.4  Control  Objective 

Hie  control  system  (i.e.  actuator  placement  and  controller 
design)  objective  is  to  minimize  tip  (node  7)  movement  in  both  the 


5  1 


x  and  y  directions.  In  this  case  the  performance  node  shapes  are 
identical  to  the  node  shapes  for  displacement  in  the  x  and  y  direction  at 
node  7.  The  objective  is  to  control  the  tip  movement  to  the  following 
mean  square  steady  state  bounds: 

.  -5  2 

x  <7x10  m 
ms  — 


-  5  2 

y  <9x10  m 
ms  — 


6.5  Open  Loop  Characteristics 

The  beam  is  modeled  as  having  a  damping  ratio  of  .001  for  each 
mode.  The  open  loop  characteristics  of  the  beam  when  subjected  to  the 
disturbances  are  listed  in  Tables  6-1  and  6-2.  The  mean  square  tip 
movement  in  the  x  direction  is  0.15  m^  and  0.048  in  the  y 
direction. 

6 .6  Controller  Design 

Actuator  placement  strategies  are  evaluated  via  the  performance 
of  a  designed  controller.  The  approach  taken  here  is  to  design  a 
constant  gain  state  feedback  regulator  and  a  constant  gain  Kalman  filter 
using  the  LQG  control  philosophy.  The  controller  is  based  on  the  primary 
modes  only.  The  primary  mode  state  equations  are 


x  *  A  +  B  u  +  ?  d 
— p  P  P~  P“ 


C  x 
P-P 


where  all  terms  are  as  defined  in  Chapter  3.  In  addition,  tip  movement 
can  be  written 
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Table  6-1  Open  Loop  Poles  and  Vibration  Frequencies 


SEAL 

-S.4106S7330E-04  5. 
-7.150722150E-04  7. 
-5.874783520E-03  5. 
-7 . 403243060E-03  7. 
-1.818135Q70E-02  1. 
-2.007872010E-02  2. 
-2.334321590E-02  2. 
-2.572767640E-02  2. 
-3.616505430E-02  3. 
-4.659422300E-02  4. 


Table  6-2  Open 


MODE 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


mis 

410684625E-01 
150716575E-01 
874780583E+00 
403239356E*00 
818134161E+01 
007871006E+01 
334320423E*01 
572766354E+01 
616503622E+01 
6594199 70E* 01 


mens 

1.000000000E-03 
1.000000000E-Q3 
1 . OOOOOOOOOE-Q3 
1.000000000E-03 
1 . 000000000E-03 
1.000000000E-03 
1. OOOOOOOQOE-Q3 
1 . OOOOOOOOOE-Q3 
1.000000000E-03 
1.000000000E-03 


epbubst 

5. 4106873 30E-01 
7.150722150E-01 
5.874783520E+00 
7.403243860E+00 
1 . 8181350 70E«0l 
2 . 007872010E*41 
2 . 334321590E+Q1 
2.572767640E+01 
3.616505430E*Q1 
4.659422300E+01 


Loop  Modal  Amplitudes 

Mean  Square  Amplitude 
55. 

44. 

15.0 

7.6 

0.59 

6.0xl0~5 

0.33 

0.016 

4.9xl0~3 

5.0xl0-4 
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where  C  is  a  matrix  containing  the  performance  node  shapes.  It  is 
L* 

assumed  that  position  sensors  are  collocated  with  the  actuators. 

The  regulator  gains  are  selected  to  minimize  the  cost  J 


where 


5  is  the  output  weighting  matrix  penalizing  tip  movement 

in  the  x  and  y  direction 


i?w  is  the  input  weighting  matrix 


p  is  a  scalar  multiplier 

The  optimal  regulator  gain  matrix,  K,  is  given  by 


*  -  p(*w)  BpS 


where  S  is  the  solution  to  the  algebraic  matrix  Riccati  equation 


0  -  -SA  -  ATS  +  pSB  (R  )  BTS  -  cJq  Ct 

p  p  K  p*>  wJ  p  L*w  I 


*n.~**~ r_ 


Note  that  each  actuator  set  will  have  a  different  B  matrix  which 

P 

requires  calculation  of  optimal  regulator  gains  for  each  set  to  be 
evaluated. 

A  constant  gain  Kalman  filter  is  implemented  as  the  estimator. 
The  constant  Kalman  gain  matrix,  G,  is  given  by 


G  =  PCT(1 

p '  n ' 


where  P  is  the  solution  of 


A  P  +•  PAT  +  TQ_ ?T  -  PCTfR  )~'c  P  =  0 

p  p  p^  n'  p 


is  the  plant  noise  variance  matrix 


R  is  the  sensor  noise  variance  matrix 
n 


Note  that  sensor  noise  is  not  included  in  the  problem  formulation.  The 

sensor  noise  variance  matrix,  R  ,  is  used  as  a  design  parameter  in  this 

n 

problem. 

The  controller  designed  for  the  SCE  beam  is  based  on  the 
following  values  for  the  weighting  matrices  and  noise  covariances. 


100  0 


800  0 


0  800 


R  - 

n 


'  '  •  .  •-  -V'-1 


p  is  a  design  variable  that  is  selected  for  each  actuator  set  such  that 
the  controller  meets  design  specifications. 


6.7  Control  System  Evaluation 

The  steady  state  performance  of  a  system  must  be  obtained  to 
evaluate  the  control  system  design.  Consider  a  system  with  the  state 
equation 


where 


x  =■  Ax  +  T  w 
—a  a— a  a — 


A  can  represent  open  or  closed  loop  dynamics 

3L 


x  is  the  appropriate  state  vector 


w  is  a  vector  of  white  noise  inputs 


T  is  the  appropriate  noise  influence  matrix 
a 


The  steady  state  covariance  matrix,  X,  is  obtained  by  solving  the 
Liapunov  equation  in  the  form^20^ 


-  _  — T  -  -  T 

a  x  +  xa  +  r  q  r 

a  a  a*n  a 


where  Q  is  the  white  noise  variance  matrix.  The  system  performance  is 
n 

defined  as  the  mean  square  value  of  the  output,  y,  given  by 


Thus,  the  mean  square  performance,  Y,  can  be  expressed  as 


Y 


-  — T 
CLXCL 


The  control  inputs  are  defined  as 


u  - 


Kx 


where  x  is  a  vector  of  primary  mode  states  for  a  full  state  feedback 
system.  For  a  system  with  an  estimator,  X  is  the  estimate  of 
the  primary  mode  states.  Thus,  the  mean  square  control  inputs,  U,  can 
be  expressed  as 

U  *  kxkt 

where  X  is  the  appropriate  state  covariance  matrix. 

Given  the  above  formulation,  the  control  system  is  evaluated  as 

follows: 

(1  )  Select  a  p  value  such  that  the  primary  mode  full  state 
feedback  system  (i.e.  regulator  but  no  observer)  meets 
specifications  for  Y. 

(2)  Calculate  U 

The  above  characteristics  can  be  compared  for  control  systems  based  on 
different  actuator  set  selections. 

6.8  Actuator  Placement 

With  the  above  preliminaries  completed,  the  actuator  selection 
procedures  developed  in  Chapters  4  and  5  can  be  applied. 
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6.8.1  Choosing  Initial  Set 

In  an  effort  to  make  this  a  non-trivial  problem,  actuators  are 
not  allowed  to  be  placed  on  the  beam  section  from  node  51  to  beam  tip,  a 
range  of  approximately  14  meters.  Since  the  beam  has  the  convenient 
property  that  the  mode  shapes  are  oriented  predominantly  in  either  the  x 
or  y  direction,  placing  actuators  in  off  axis  directions  can  be  ruled 
out.  The  actuators  considered  are  of  the  point  force  type.  The  initial 
set  consists  of  x  and  y  direction  actuators  located  at  nodes  26  through 
51.  This  gives  a  total  of  52  possible  actuator  locations.  Actuator 
numbers,  location,  and  directions  of  influence  are  tabulated  in  Appendix 
A. 


6.8.2  Primary  Mode  Selection 

The  first  stop  n  the  reduced  order  controller  design  is  to  select 
the  primary  modes.  For  this  system,  the  modes  with  the  most  contribution 
to  tip  movement  will  be  selected.  The  overall  tip  deflection,  m,  is 
given  by 


m 


n 


i  X, 


TIP 


+  y 


2 

TIP 


The  approximate  contribution  of  each  mode  to  the  tip  deflection  can  be 
calculated.  These  values  are  listed  in  Table  6-3  for  the  uncontrolled, 
disturbed  beam.  Based  on  these  values,  modes  7,  8,  9,  and  10  are 
selected  as  primary  modes. 
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Table  6-3.  Mean  Square  Contributions  To  Tip  Deflection 


Mode 

x  (m2) 
ms  ' 

y  f  m2} 
ms  ^  ' 

m  f  m2 ) 
ms 1  1 

7 

0.150 

1  . 1  xl  0_1  7 

0.150 

8 

8.0X10*1 7 

0.0459 

0.0459 

9 

2.2x10~14 

2 .23x1 0~3 

2 .23x1  0_3 

10 

1  .1  3x1  0-3 

2.5x10~14 

1  .1  3x10~3 

1 1 

1  .1x10_1  3 

2.8x1 0~3 

2.8x10-5 

12 

3.2x1 0~23 

2.9x10_1 3 

2.9x1  0_1 3 

1  3 

1 .5x1 0~3 

1  .1x10~12 

1 .6x1 0-3 

14 

1  .1x10-12 

3 .5x1 0~8 

3.5x1 0~8 

15 

3 .2x1  0_1  5 

9.0x10~8 

9.0x1 0~8 

16 

9.6x1 0~9 

1 .8x1 0*1 6 

9.6x1 0~9 

As  noted  in  Chapter  4,  the  selection  of  the  elements  in  the  least 
squares  weighting  matrix,  5,  is  based  on  the  contribution  of  each  primary 
mode  to  overall  performance.  Thus,  the  elements  of  the  D  matrix  are 
obtained  from  Table  6.3.  The  D  matrix  obtained  is 

D  »  diag(0.150,  0.0459,  2.23x10~3,  1.13x10-3,  0,  0,  0,  0,  0,  o) 


6.8.3  Placement  via  Maximum  Norm 

Using  the  maximum  norm  selection  procedures  developed  in  Chapter 
5,  the  actuator  set 


51  -  36  -  52  -  33 
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is  obtained.  The  R  matrix  obtained  from  the  QR  factorization  procedure 


5  l 

36 

5  2 

33 

1) 

1.57398E-02 

-1.10470E-10 

-3. 56998E-10 

5.838846-03 

2) 

0.0 

7.43986E-03 

2.99312E-03 

2.57983E-09 

3 1 

0.0 

0.0 

4.56074E-03 

-1.94172E-09 

4) 

0.0 

0.0 

0.0 

4. 33872E-03 

The  appropriate  oi  vectors  for  the  x  and  y  directions  are 


—  - 

-1 

-7 

-1  .18x10 

-2.73x10 

-7 

-1 

2.13x10 

a  s 

9.42x10 

3.06x1 0“B 

-y 

-3.20x10“' 

-1 

-7 

8.61x10 
■» .  — 

8.38x10 

Actuator  locations  are  plotted  on  Figure  6-2. 

6.8.4  Placement  via  Residual  Matching 

In  this  section  the  Residual  matching  procedure  with  no  spillover 
considerations  (i.e.  4*1.0)  is  used  to  place  actuators.  The  selected  set 
for  each  direction  is 


X 

Z_ 

51 

52 

39 

34 

34 

51 

50 

39 

This  set  of  5  actuators  (note  duplications)  is  reduced  to  4  actuators  by 
the  maximum  norm  procedure.  The  resulting  final  set  is 
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51  -  34  -  39  -  52 


for  which  the  R  matrix  is 


51 

34 

5  2 

3  9 

1) 

1.57398E-02 

-1.10470E-i0 

-3.5699SE-10 

8.91708E-03 

2) 

0.0 

7.439S6E-03 

2.99312E-03 

1.27244E-09 

3) 

0.0 

0.0 

4.56074E-03 

-1.84773E-09 

4) 

0.0 

0.0 

0.0 

3. 93113E-03 

and  the  a  vectors  are 


-1 

-7 

-3.37x10 

-4.87x10 

-7 

-1 

3.42x10 

a  * 

9.42x10 

4.91x1 0~° 

-y 

-3.20x10 

-1 

-7 

9.51x10 

9.24x10 

Actuator  locations  are  plotted  on  Figure  6-3. 


6.8.5  Comparing  the  Two  Methods 

The  above  actuator  selections  are  based  on  finding  the  actuator 
set  with  the  best  control  authority  over  the  system.  Each  set  matches 
the  intuitive  notion  that  actuators  should  be  placed  near  peaks  in 
primary  mode  shapes  to  give  good  control  authority.  A  full  state 
feedback  controller  based  on  the  primary  modes  is  used  to  compare  the  two 
actuator  sets. 

For  the  maximum  norm  set  (51-36-52-33),  the  design  parameter  p  is 
chosen  as  .0001.  this  gives  a  mean  square  performance  of 

x  -  6.91  x  10~5  m2 

ms 
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-5  2 

y  =  8.69x10  m 

ms 


The  control  input  covariance  matrix  is 


1) 

2) 

31 

4) 


51 

3.520559E+01 

2.915335E«-01 

-1.166614E-C5 

-X.747576E-06 


33 

2.9X5335E+01 

4.683268E*01 

-1.122090E-05 

2.075952E-0S 


52 

-1.166614E-05 

-1.122090E-05 

1.136125E+02 

1.848S57E+01 


36 

-X.747576E-06 
2.075952E-0S 
1.848557E+01 
8 . 968870E+01 


A  scalar  measure  of  the  total  force  input  is  the  trace  of  the  covariance 
matrix,  Tr{u}.  For  this  set 


Tr{5}  =  285.33  nt2 

For  the  residual  matching  set  (51-34-52-39),  the  design  parameter 
p  is  again  chosen  as  .0001 .  This  gives  a  mean  square  performance  of 


-5  2 

x  ■  6.71  x  10  m 

ms 


-5  2 

y  *  8.64  x  10  m 

ms 


The  control  input  covariance  matrix  is 


51 

39  . 

.52  . 

34 

1) 

2.943699E+01 

3. 093736E+01 

-8.171109E-06 

X.007317E-0S 

2) 

3.093736E+01 

4.944458E+01 

-8.596048E-06 

3.242598E-05 

3) 

-6.171109E-06 

-8.S96048E-06 

1.087050E+02 

1.167082E«01 

4) 

X.  00731 7E-05 

3.242S98E-05 

1.167082E+Q1 

9.478019E*01 
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For  which 


Tr{u}  =  282.37  nt2 

The  actuator  set  selected  by  intuitive  methods  (i.e.  peaks  of 
mode  shapes)  is  (5 1-52-37-38) .  Selecting  the  design  parameter  p  as  .0001 
gives  a  mean  square  performance  of 


x  ■  6.68x10 

ms 


y  -  8.82  x  10 

ms 


The  control  input  covariance  matrix  is 


51 

52 

3  7 

38 

1) 

3.102040E+01 

-8.179020E-06 

2.995857E+01 

1.631726E-C6 

2) 

-8.179020E-06 

1 . 171245E+02 

-7.931640E-06 

2.891552E+01 

31 

2. 995G57E+01 

-7.931640E-C6 

4.82691 9E+01 

2.S66040E-C5 

4) 

1.631726E-06 

2.S91552E+01 

1.866040E-05 

8.663411E+01 

for  which 

Tr{u}  -  282.65 


The  above  results  show  the  residual  matching  procedure  yields  a 
better  actuator  set  than  the  maximum  norm  procedure  in  the  sense  that 
slightly  better  performance  was  obtained  with  smaller  control  force 
inputs. 
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6.8.5  Placement  via  Spillover  Considerations 

The  residual  matching  procedure  was  developed  with  the  capability 
of  placing  actuators  to  reduce  spillover  to  selected  secondary  modes. 

The  selection  of  secondary  modes  is  not  a  clear  cut  procedure.  For  this 
study,  residual  modes  that  appear  most  sensitive  to  spillover  effect  will 
be  designated  as  secondary  modes. 

To  determine  those  modes  sensitive  to  spillover  effect,  the 
actuator  set  51-39-52-34  is  used  for  a  reduced  order  controller  design 
(regulator  and  estimator)  that  is  applied  to  the  full  10  mode  system. 

With  the  design  parameter  p  =  .0001  as  utilized  above,  the  system  is 
unstable.  Pole  locations  are  listed  in  Table  6-4.  Mode  13  is  the 
unstable  mode,  and  the  pole  for  mode  11  has  been  shifted  very  close  to 
the  imaginary  axis.  Designating  modes  11  and  1 3  as  secondary  modes,  an 
attempt  is  made  to  select  a  set  of  actuators  that  will  meet  system 
specifications  while  maintaining  a  stable  system. 

Using  the  residual  matching  procedure  with  a  limiting  factor 
i” 0.9,  the  actuator  set  51-35-52-36  is  obtained.  With  a  design  parameter 
of  p”9.0x10-5  in  the  full  state  feedback  system,  the  mean  square 
performance  is 


x  *  6.39  x  1 0 

ms 


y  ■  8.19x10 

ms 


The  control  input  covariance  matrix  is 


5  1 

35 

52 

36 

1) 

3.432231E+01 

3.oeeis6E*oi 

-9.140142E-06 

4.674397E-06 

2) 

3.0«ai«6E*01 

5.001754E+01 

-7.502096E-06 

2.973809E-05 

3) 

-9.140142E-06 

-7.502096E-06 

1.186631E+02 

1.915437E*01 

4) 

4.674397E-06 

2.973809E-05 

1.91S437E+01 

9.497985E+0X 

67 


for  which  Tr{u}  =*  298.18.  This  indicates  slightly  higher  force  inputs 
are  required  to  match  previous  performance  levels.  However,  when  this 
set  is  used  in  the  design  of  the  reduced  order  controller,  the  full  order 
controlled  system  remains  stable.  Thus,  in  this  simple  example,  a 
reasonable  selection  of  actuator  locations  has  been  obtained  which 
achieve  the  control  objectives  despite  the  spillover  effect. 


CHAPTER  7 


SUMMARY  AND  CONCLUSIONS 


7.1  Summary 

The  purpose  of  this  study  was  to  develop  a  method  of  actuator 
placement  to  allow  vibration  control  of  large  spce  structures. 

Previously  developed  techniques  have  tended  to  be  computationally 
intensive  and  thus  are  not  easily  applicable  to  high  order  systems.  The 
methods  for  actuator  placement  presented  in  this  study  are  readily 
extended  to  high  order  sysems  and  shewed  much  promise  when  applied  to  a 
simple  large  space  structure. 

With  the  system  equations  written  in  modal  coordinates,  vectors 
of  displacements  associated  with  both  actuator  and  system  performance 
could  be  developed.  These  were  termed  node  shapes.  The  actuator 
placement  problem  was  formulated  as  a  least  squares  problem  where  a  set 
of  actuator  node  shapes  from  an  initial  set  of  n  actuators. 

The  first  technique  initially  selected  the  actuator  node  shape 
with  the  maximum  norm.  Subsequent  actuator  node  shapes  were  chosen  based 
on  the  component  of  the  node  shape  orthogonal  to  the  space  spanned  by  the 
previously  selected  set.  The  orthogonal  conqponent  of  each  node  shape  was 
calculated  and  the  component  with  the  maximum  norm  was  added  to  the 
actuator  set.  This  method  was  termed  the  maximum  norm  technique. 

The  second  technique  initially  selected  the  actuator  node  shape 
with  the  largest  absolute  projection  along  the  performance  node  shape. 

The  least  squares  residual  vector  could  then  be  calculated  and  subsequent 
selections  were  based  on  the  actuator  node  shapes  absolute  projection 
along  the  residual  vector.  This  method  was  termed  the  residual  matching 
technique. 


Application  of  these  techniques  to  a  simple  beam  structure 
demonstrated  that  actuator  locations  chosen  were  near  peaks  of  primary 
mode  shapes,  which  is  the  intuitive  concept  for  actuator  location.  The 
maximum  norm  technique  appeared  to  choose  actuator  locations  that  were 
more  orthogonal  than  those  attained  by  intuitive  placement.  The  residual 
matching  technique  appeared  to  choose  actuator  locations  that  were  more 
advantageous  to  applying  control  inputs  to  all  the  primary  modes  with  a 
single  actuator,  rather  than  an  individual  mode  as  in  intuitive 
placement.  In  other  words,  the  residual  matching  technique  sacrificed 
some  actuator  orthogonality  in  order  to  attain  slightly  better  control 
authority  over  the  primary  modes.  When  the  control  forces  required  to 
meet  a  specific  performance  objected  were  examined,  the  residual  matching 
actuator  set  gave  approximately  the  same  performance  as  the  intuitive 
actuator  set,  and  both  performed  better  than  the  maximum  norm  actuator 
set. 

Position  sensors  were  assumed  to  be  colocated  with  the  actuators 
and  a  Kalman  filter  was  added  to  the  control  design.  Hie  actuator 
locations  used  were  those  selected  by  the  residual  matching  technique. 
Analysis  of  the  system  poles  shewed  that  one  of  the  poles  was  shifted 
into  the  right  half  s-plane  by  spillover  effect.  Hie  residual  matching 
technique  was  modified  to  take  spillover  effect  into  account  in  the 
selection  process  and  actuator  locations  were  reselected.  Hie  resulting 
set  required  higher  control  force  levels  to  meet  performance 
specifications  in  the  full  state  feedback  system,  but  when  the  system 
poles  were  calculated  with  the  Kalman  filter  in  place,  all  poles  remained 
in  the  left  hand  s-plane. 

7.2  Conclusions  and  Recommendations 

Hie  above  application  demonstrates  that  both  selection  techniques 
show  promise  as  actuator  placement  techniques.  In  this  case  the  residual 
matching  technique  appears  to  be  superior  since  the  actuator  set  chosen 
by  this  technique  allows  a  full  state  feedback  controller  to  operate  with 


lower  input  force  levels  than  the  set  selected  by  the  maximum  norm 
method.  In  addition,  the  residual  matching  technique  has  been  modified 
to  select  actuators  so  as  to  limit  spillover  into  selected  modes.  The 
application  problem  demonstrates  that  proper  actuator  placement  can  lead 
to  a  stable  system  that  is  capable  of  meeting  performance  specifications. 

the  example  structure  is  geometrically  simple  and  is  represented 
by  a  low  order  model.  Application  of  the  selection  techniques  to 
structures  with  more  complex  geometry  and  models  of  higher  order  is 
required  to  verify  the  suitability  of  these  placement  techniques. 


APPENDIX  A 
SCE  BEAM  DATA 


Table  A-l  Truss  Characteristics 


Truss  characteristics 

El  pitch 

2.0  X  107  N-m2 

El  roll 

1.3  X  107  N-m2 

Length 

100m 

i 

2.0m 

h 

2.83m 

Table  A-2  Node  Locations 


NODE 

X(m) 

Y(m) 

Z(m) 

1 

0.0 

0.0 

0.0 

2 

0.0 

0.0 

0.0 

3 

0.0 

-0.173740005 

£♦01 

0.0 

4 

-0 . 347679996 

E+01 

-0.110235977 

Em 

0.0 

5 

0.0 

0.173740005 

£♦01 

0.0 

6 

-0.347679996 

£♦01 

0.110235977 

£+01 

0.0 

7 

0.0 

0.0 

0.100457199 

E+03 

a 

0.0 

-0.110235977 

E+01 

0.0 

9 

0.0 

-0.110235977 

E+01 

0.457199991 

E+00 

10 

0.0 

-0.110235977 

E+01 

0.245718956 

E+01 

11 

0.0 

-0.110235977 

E+01 

0 . 3484SS045 

E+01 

12 

0.0 

-0.110235977 

E  +  01 

0.445718956 

E+01 

13 

0.0 

-0.110235977 

E  +  01 

0.579119968 

E+01 

14 

0.0 

-0.110235977 

E+01 

0.645718956 

E+01 

15 

0.0 

-0.110235977 

E+01 

0.845718002 

E+01 

16 

0.0 

-0.110235977 

E+01 

0.929487038 

E+01 

17 

0.0 

0.110235977 

E+01 

0.0 

la 

0.0 

0.110235977 

E+01 

0.457199991 

E+00 

19 

0.0 

0.110235977 

E+01 

0.245718956 

E+01 

20 

0.0 

0.110235977 

E+01 

0.348488045 

E+01 

21 

0.0 

0.110235977 

E+01 

0.445718956 

E+01 

22 

0.0 

0.110235977 

E+01 

0.579119968 

E+01 

23 

0.0 

0.110235977 

E+01 

0.645718956 

E+01 

24 

0.0 

0.110235977 

E+01 

0.845718002 

E+01 

25 

0.0 

0.110235977 

E+01 

0.929487033 

E+01 

26 

0.0 

0.0 

0.457199991 

E+00 

27 

0.0 

0.0 

0.245720005 

E+01 

28 

0.0 

0.0 

0.445720005 

E+01 

29 

0.0 

0.0 

0.645720005 

E+01 

30 

0.0 

0.0 

0.845720005 

E+01 

31 

0.0 

0.0 

0.104572001 

E+02 

32 

0.0 

0.0 

0.124572001 

E+02 

33 

0.0 

0.0 

0.144572001 

E+02 

34 

0.0 

0.0 

0.184571991 

E+02 

35 

0.0 

0.0 

0.224571991 

E+02 

36 

0.0 

0.0 

0.264571991 

E+02 

37 

0.0 

0.0 

0.304571991 

E+02 

3a 

0.0 

0.0 

0.344571991 

E+02 

39 

0.0 

0.0 

0.384571991 

E+02 

40 

0.0 

0.0 

0.424571991 

E+02 

41 

0.0 

0.0 

0.464571991 

E+02 

42 

0.0 

0.0 

0.504571991 

E+02 

43 

0.0 

0.0 

0.544571991 

E+02 

44 

0.0 

0.0 

0.584571991 

E+02 

45 

0.0 

0.0 

0.624571991 

E+02 

46 

0.0 

0.0 

0.664571991 

E+02 

47 

0.0 

0.0 

0.704571991 

E+02 

48 

0.0 

0.0 

0.744571991 

E+02 

49 

0.0 

0.0 

0.784571991 

E+02 

50 

0.0 

0.0 

0.824571991 

E+02 

51 

0.0 

0.0 

0.664571991 

E+02 

52 

0.0 

0.0 

0.904571991 

E+02 

53 

0.0 

0.0 

0.944571991 

E+02 

54 

0.0 

0.0 

0.984571991 

E+02 
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Table  A- 3  Element  Connection  Data 


element 

10001 

10002 

10003 

10004 

10005 

10006 

10007 

10008 
20001 
20002 

20003 

20004 

20005 

20006 
20007 
20003 

30001 

30002 

30003 

30004 

30005 

30006 

30007 
30003 

30009 

30010 

30011 

30012 

30013 

30014 

30015 

30016 

30017 

30018 

30019 

30020 

30021 

30022 

30023 

30024 

30025 

30026 

30027 
30023 
30029 

111 

112 

121 

211 

212 

221 

31011 

31012 

31013 

31014 

31015 

31021 

31022 

31023 

31024 

31025 


NODEl 

NODE2 

8 

9 

9 

10 

10 

11 

11 

12 

12 

13 

13 

14 

14 

15 

15 

16 

17 

13 

18 

19 

19 

20 

20 

21 

21 

22 

22 

23 

23 

24 

24 

25 

26 

27 

27 

28 

28 

29 

29 

30 

30 

31 

31 

32 

32 

33 

33 

34 

34 

35 

35 

36 

36 

37 

37 

38 

33 

39 

39 

40 

40 

41 

41 

42 

42 

43 

43 

44 

44 

45 

45 

46 

46 

47 

47 

48 

48 

49 

49 

50 

5Q 

51 

51 

52 

52 

53 

53 

54 

54 

7 

3 

10 

3 

12 

4 

10 

5 

19 

5 

21 

6 

19 

26 

9 

27 

10 

23 

12 

29 

14 

30 

15 

26 

18 

27 

19 

23 

21 

29 

23 

33 

24 

75 


Table  A-4 

Actuator 

Locations 

ACTUATCS 

BISECTION 

NODE 

1) 

TX 

26 

2) 

7Y 

26 

3) 

TX 

27 

4) 

TY 

27 

5) 

TX 

28 

6) 

TY 

28 

7 ) 

TX 

29 

8) 

TY 

29 

9) 

TX 

30 

10) 

TY 

30 

11) 

TX 

31 

12) 

TY 

31 

13) 

TX 

32 

14) 

TY 

32 

15) 

TX 

33 

161 

TY 

33 

17) 

TX 

34 

IS) 

TY 

34 

19) 

TX 

35 

20) 

TY 

35 

21) 

TX 

36 

22) 

TY 

36 

23) 

TX 

37 

24) 

TY 

37 

25) 

TX 

38 

26) 

TY 

38 

27) 

TX 

39 

2S) 

TY 

39 

29) 

TX 

40 

30) 

TY 

40 

31) 

TX 

41 

32) 

TY 

41 

33) 

TX 

42 

34) 

TY 

42 

35) 

TX 

43 

36) 

TY 

43 

37) 

TX 

44 

36) 

TY 

44 

39) 

TX 

45 

40) 

TY 

45 

41) 

TX 

46 

42) 

TY 

46 

43) 

TX 

47 

44) 

TY 

47 

45) 

TX 

48 

46) 

TY 

48 

47) 

TX 

49 

43) 

TY 

49 

49) 

TX 

50 

50) 

TY 

50 

51) 

TX 

51 

52) 

TY 

51 

53) 

TX 

52 

54) 

TY 

52 

55) 

TX 

53 

56) 

TY 

53 

57) 

TX 

54 

53) 

TY 

54 
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Figure  A- 2  Mode  8  Mode  Shape 
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Figure  A- 5  Mode  11  Mode  Shape 
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Figure  A-9  Mode  15  Mode  Shape 
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